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Abstract. We study time-reversal symmetry in dynamical systems with finite phase space, 
with applications to birational maps reduced over finite fields. For a polynomial automor- 
phism with a single family of reversing symmetries, a universal (i.e., map-independent) distri- 
bution function 1Z(x) = 1 — e~ x (l + x) has been conjectured to exist, for the normalized cycle 
lengths of the reduced map in the large field limit [2D]. We show that these statistics corre- 
spond to those of a composition of two random involutions, having an appropriate number 
of fixed points. This model also explains the experimental observation that, asymptotically, 
almost all cycles are symmetrical, and that the probability of occurrence of repeated periods 
is governed by a Poisson law. 



1. Introduction 

This paper is concerned with time-reversal symmetry in dynamical systems with finite 
phase space, with applications to rational maps over finite fields. The concept of reversibility 
originated in the theory of smooth maps and flows on manifolds El H]. In this paper, a 
map L is said to be reversible if it is the composition of two involutions G and H = L o G. 
The involution G conjugates the map to its inverse, namely 

(1) GoLoG- 1 = L~ l . 

(For background information on reversibility and its generalization, see [T^l 118] , and 
references therein.) The nature of the fixed sets of G and H, respectively FixG and FixH, 
plays a key dynamical role in organizing the dynamics of a reversible map. 

The reversibility property can be interpreted in purely algebraic terms, and in particular it 
applies to reversible algebraic maps L of any dimension n over any field K, namely reversible 
maps defined by rational functions with coefficients in K . If the field K = ¥ q is a finite field 
with q elements — where q is the power of a prime number p — we are led to the study of 
reversibility on a phase space with N = q n points. (For background reference on finite fields, 
see [13].) More meaningful is the case in which K is an algebraic number field, e.g., K = Q. 
The same map L may then be reduced over infinitely many distinct finite fields, of increasing 
size. 

We are concerned here with the case where the reductions of the component involutions G 
and H are involutory permutations of the phase space of N = q n points, so that the reduction 
of their composition L yields a permutation of this space. In reduction, the fixed sets Fix G 
and Fix H will now be finite. Examples arising from dimension n = 2 are the reversible planar 
polynomial automorphisms [H [8] which include the well-known area-preserving Henon map 
Li = ffioGi, 



(2) G\ : (x,y) i— > (y, x) FL\ : (x, y) i-> (x, -y + x 2 + a) a £ K. 
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If K = F p , we find 

#FixG x = #FixHi = p. 

Reversible polynomial automorphisms in any dimension provide other examples. But cer- 
tain rational maps in n dimensions can induce permutations, e.g., L 2 = H 2 o G 2 of [T7] 

G<z:(x,y,z) 1— ► (x + e(2y — k)(z + e(y — k)), k — y, z + e(2y — k)) e,k&K 
H 2 :(x,y,z) ^ (y (2 - 2x/F + x 2 /F 2 ), x/F, -z) F = l + (l-y) 2 . 
If K = F p , with p = 3 (mod 4), the denominator in H 2 is non- vanishing and 

#FixG 2 =p 2 , #FixH 2 =p. 

The representability of G and H over a finite field ¥ q imposes some restrictions on the 
parameters present in the involutions. For instance, if in ([2]), we take a G Q, say, we are led 
to consider the set P a of the prime divisors of the denominator of a. Then the involution H\ is 
representable over any finite field whose characteristic is not in P a , such as the field of integers 
modulo a prime p P a . Additionally, one must verify that G and H remain involutions when 
reduced to a finite field. This is an instance of the property of good reduction (for a discussion 
of reduction in arithmetic dynamics, see \21\ chapter 2]). 

Representability assumed, the reduction of an algebraic map can be associated with in- 
teresting asymptotic phenomena, concerning the statistical behaviour of the periodic orbits 
in the large field limit. A main theme in this area of research is to ascertain whether such 
phenomena are universal, namely independent from the map, depending only on its struc- 
tural properties e.g., reversibility. Furthermore, in reduction, the original dimension n of the 
underlying algebraic map ceases to be important as the finite phase space does not inherit 
the topology of the continuum M n or C n . 

If the reduction of a reversible algebraic map is a permutation, it has been shown [20J that 
over finite fields, reversibility manifests itself combinatorially. In particular, the number of 
symmetric periodic orbits — those invariant under G — equals 

(3) (#FixG + #FixH)/2. 

Strong experimental evidence [201 M suggests that the symmetric periodic orbits dominate 
the statistics over the asymmetric onesj, so that (J3J) asymptotically counts the number of 
periodic orbits in the permutation. Furthermore, numerical experiments [20 \ [9] suggest there 
exists an asymptotic (large fields) distribution of the periods of the orbits, given by 

(4) K(x) = 1 - e~ x (l + x) x>0. 

(In these experiments, 'large fields' refers to the limit ¥ p , p — > 00; there are other asymptotic 
regimes, see remark (3) in section SJ) The distribution dH) describes the limiting probability 
of the set of points belonging to cycles with scaled period not exceeding x. This distribution 
is believed to be universal, within the class of reversible maps with a single time-reversal 
symmetry. By contrast, maps which are not reversible appear to behave like random permu- 
tations [20]; the associated distribution, also conjectured to be universal, is markedly different 
from that of reversible maps. Finally, there are asymptotic period distributions for inte- 
grable systems [HI [10]. Even though these distributions are map-dependent, they all feature 
a distinctive 'quantization' of periods. With appropriate scaling, the allowed periods occur 

^This is in sharp contrast with the case of maps with real or complex coordinates [14] . 
2 The inadequacy of random permutations as a model for reversible maps on a discrete space was noted as 
far back as [To] . 
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at the reciprocals of the natural numbers, resulting in step-like distribution functions. These 
various asymptotic phenomena have led to the development of simple and effective tests for 
detecting integrability and reversibility in algebraic mappings [19\ [20l [TO] . 

In this paper we propose a combinatorial model for the orbits, over a finite field, of a map 
with a single family of reversing symmetries. Our model consists of the composition of two 
random involutions. The size N of the space and the numbers j^FixG and #FixH of fixed 
points of each involution are the parameters of the model. By letting N increase, we build a 
sequence of probability spaces. Under very general conditions, we derive the distribution (@J 
found in [201 [9]. 

We prove the following theorem, which appears in section [2] as theorems Q] and El 

Theorem A. Let (G, H) be a pair of random involutions of a set O with N points, and let 
g = j^FixG and h = j^FixH. Let TZn(x) be the expectation value of the portion of Cl 
occupied by cycles of H oG with period less than 2xN/ (g + h), computed with respect to the 
uniform probability If, with increasing N, g and h satisfy the conditions 

(5) lim g(N)+h(N) = oo lim 9{N) + h(N) =Q 

then for all x > 0, we have the limit 7Zn(x) — > TZ(x), the distribution function in 
Moreover, almost all points in f2 belong to symmetric cycles. 

This result is formulated in terms of the sequence g + h, and not of the two sequences 
individually. The conditions on the growth rate of g+h are quite mild, hence the wide scope of 
applicability of the theorem. If g and h do not grow at the same rate, an interesting dynamical 
phenomenon occurs: asymptotically, almost all periodic orbits have even period (corollary [9]) . 
For n-dimensional maps over the finite field ¥ q , one has N = q n , while g and h will typically 
grow algebraically (see remark 1 of section [4]) . As mentioned above, increasing the dimension 
does not introduce new difficulties, because the absence of a significant topological structure 
makes the combinatorial model effectively dimension- independent. All that is needed is to 
keep track of the possible cardinalities of the fixed sets of the two involutions. 

In [20], experimental evidence from reductions of the Henon map ([2]) suggested that the 
probability of occurrence of several cycles with the same period followed a Poisson law, with 
parameter depending exponentially on the (rescaled) period. These findings motivate the 
study of the statistics of the occurrence of repeated periods in the random involutions model. 
In section [3] we consider this question, and prove the following result. 

Theorem B. With the notation of theorem A, let fi(t, i) be the probability that (G, H) has 
i cycles of period t. Assume that the sequence 



(6) /(AO = <i 

has the (possibly inhnite ) limit 



if t is odd 

N 

9 2 + h 2 

it t is even 

2N 



c= lim f(N) 

N^oo 

and let 

p) y.iizmm-w 
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Then, if c ^ we have, as N — > oo 

, —n,OL % Ic if t is constant 

(8) MM) ~ e t «= S _ vftl ... 

I e ' if y is constant 

while if c = we Lave [i(t, i) ~ <5j, where 5 is Kronecker's delta. 

As it will appear in the proof of theorem B, in the case of constant y, for convergence to a 
Poisson distribution it is only required that f(N) /> 0, not the actual existence of a non-zero 
limit c. We conclude section 3 by showing the agreement between the prediction of theorem 
B with the actual repetitions i) for a small t in the Henon map. The latter are governed 
by the roots (mod p) of certain polynomial functions of one variable, with a connection to 
the Galois groups of these polynomials. 

In section 2] we discuss the question of the asymptotic (large field) growth rate of the 
parameters associated with a rational map over a finite field, namely the numbers g and h 
of fixed points of the involutions. In this case N is the power of a prime number. We show 
that the growth conditions of theorem A on g + h are easily satisfied, while f(N) in theorem 
B depends algebraically on N, and all scenarios considered in the theorem can be realized. 
Finally, we also sketch how the random involution model might be extended to the case of 
reversible rational maps with singularities. 



2. Composition of random involutions 

Let N be a positive integer, and let g, h be integers in the range 1 < g, h < N such that 
N — g and N — h are both even. We consider ordered pairs (G, H) of random involutions on 
a set of N points, which fix g and h points, respectively. Thus g = j^Fix G, h = j^Fix H, 
and the cycle decomposition of G consists of g fixed points and (N — g)/2 two-cycles, and 
similarly for H. For each pair (G, H ) of involutions we consider the composition L = H o G, 
which is a reversible permutation of f2. 

Let G and H be the sets of all involutions with the given parameters. We regard the space 

(9) E = E(g,h, N) = G x H 

as a probability space, with the uniform probability. We define the random variable Ft : E — > 
R to be the fraction of the space O occupied by the i-cycles of the map H o G, namely 

(10) Pt = — #{x G 0, : x has minimal period t under H o G} t = 1, 2, 

For given N, the sequence (Pt) has only a finite number of non-zero terms. 
The main object under study is the distribution function 

\xz\ 

(11) n N (x) = J2( p t) x>o 

t=i 

where [ - J is the floor (greatest integer part) function, z is a scaling parameter to be specified 
below, and the average (•) is computed with respect to the uniform probability on E (i.e., 
(Pt) = (#E) _1 J2e Pt)- If \xz\ = 0, the sum is empty, and TZn(x) is defined to be zero. 

The asymptotic parameter is N. We shall assume that g = g(N) and h = h(N) are positive 
integer sequences satisfying ([5]) . These sequences determine the sequence of probability spaces 
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to be considered. To get non trivial asymptotics, the scaling parameter z must grow with N 
at an appropriate rate. To this end, we define the rational sequence 

(12) ^ = WTTHN) 

which, due to ([5]), diverges to infinity. The quantity z(N) is the ratio of the size of the 
phase space to the number of symmetric cycles — cf. equation (|3|). Its full significance will be 
clarified at the end of this section (corollary [8]) . 

The main purpose of this section is to prove theorem A of the introduction, which is an 
amalgamation of theorems Q] and [7J 

Theorem 1. The distribution function TZ]\r defined in with scaling sequence z(N) given 
by ( f71|) . admits the following limit 

lim K N {x) = K{x) = 1 - e~ x (l + x) x>0. 

N— >oo 

It turns out that to establish theorem Q] it suffices to consider symmetric orbits only, 
because they provide the dominant contribution to the probability (Pt) in equation (fTTI) - 
see theorem [7J below. Accordingly, we separate out the contribution deriving from symmetric 
and asymmetric cycles as follows 

(13) P t = P t (s) + P t {a) . 

We begin with symmetric cycles of odd period t = 2k — 1. It is known [U [5j [TT1 [14], [18] 
that a cycle is of this type if and only if it corresponds to an orbit of L starting from Fix G 
and reaching Fix H for the first time after k iterations. Here and below, we follow the 
cycle dynamics keeping track of the action of G and H separately, which leads to a maximal 
sequence of 2/c points 

(14) (xi,yi, . . . ,x k ,y k ) yj = G(xj), x j+ i = H(y j ) k > 1, j = l,...,k 

with the property that all Xj are distinct. We call such a sequence a k-arc. In the case of odd 
period, we have the condition 

(15) xxeFixG, y k GFixH 

(hence y\ = x\) together with the requirement that the sequence contains t = 2k — 1 distinct 
points (see figured]). We must now determine the total number of fc-arcs of this type generated 
by the elements of E. 

Throughout the paper, the underlined exponent will denote the falling factorial powers 

(16) n fl =l rt k = n{n- l)(n- 2) ••• (n - o + 1), a > 1. 
The information we require is contained in the following two lemmas. 
Lemma 2. Let E = E(g, h, N) be as in equation (0). Then 

1 (AH) 2 



(17) #E(9,h,N) 



2 N~( g +h)/2 g \ h \ ({N - g)/2)\{{N -h)/2)V 



Proof. The cardinality of the set of all involutions on N points, with g fixed points, is given 
by 

'N\ (N-g)\ 



9 J 2W-e)l*{(N-g)/2)\ 
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FIGURE 1. Construction of a 4-arc corresponding to a 7-cycle. The solid and dashed 
lines denote the action of G and H , respectively. The orbit of L = H o G from Fix G 
to Fix H is represented by white circles; the black circles represent the return orbit 
from Fix H to Fix G. 



where the first term gives the number of ways of choosing the fixed points, and the sec- 
ond counts the number of ways of partitioning the remaining points into unordered pairs. 
Multiplying the expression above with the corresponding expression for h fixed points, we 
obtain 

1 (AH) 2 



as desired. 



#E{g,h,N) 



2 N~( 9 +h)/2 g \ h \ _ g y 2 )\ ((JV - h)/2)\ 



□ 



(s) 

Lemma 3. Let _P t be as above. The following holds 
(18) 



/ P W v. 2k ~ 1 



where 



(19) 



N 



2fc-l 



#E(g- l,h- l,N-2k + l) 



#E(g,h,N) 



(gh 

N 2 



gh 



k-2 



k > 1 



ifk = l 



Yl(N-g-2j)(N-h-2j) ifk>l. 



j=0 



Proof. Let (xi,yi,... ,Xk,yk) be the required fe-arc (see equation p^|) and figure [T]). The 
initial point x\ belongs to FixG, so that y\ = x\. The first half of the cycle consists of the 
points xi, . . . ,Xk, the second half of the points y^, . . . , y2, giving 2k — 1 distinct points in total. 
Hence the desired number of arcs is 



2k-l 



N- 

Each A;-arc occurs with multiplicity determined by the unconstrained action of the pairs of 
involutions on the rest of the space, which is given by 



#E(g-l,h-l,N-2k + l). 
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Thus the average value of -F^fe-i ^ s gi ven by 

(s) 2fe-l 2fc _! #E(g - 1, ft - 1, jV - 2fc + 1) 
\^2k-i) N #E(g,h,N) 

Finally, the product formula for <?2fc-i is obtained from lemma [2] with a straightforward 
computation. d 

By analogy with equation (jlip , we define the distribution function 1Z^'°^ for odd symmetric 
cycles as 

L(xz+1)/2J 

(20) ^* 0) (x)= £ (P«_!) ^>0. 

fc=l 

We have 

Theorem 4. For any sequences g, h satisfying |5p> z as in flfy) . we have, as N —* oo 

,( s ,o)^ _ 2gft 
(5 + /i) S 

Note that the limit distribution is symmetrical in g and /i. To prove theorem 0] we need a 
lemma. 



Lemma 5. Let g(N),h(N) be as in |5]), let m = m(N) be a positive integer sequence, and 
let 

J fc=i 

Then, as N — > oo 

(21) 5 m ~-|[l-e-" MS (m« + l)] « = 

Proof. For fixed N, m and A are fixed, and the quantity S m is the sum of an arithmetico- 
geometric progression, with value 

S m = (fTA)2 1 1 + A " A m [2m(l - A) + 1 + A] }. 

As N — > oo, we have, from ([5]) 

A ~ 1 (1 — A) ~ k m (A) ~ — k 

so we find 

S m ~ \ [2-e mla ^(2mK + 2)] ~ A [i _ e -™( m « + 1)] 
as desired. D 



Proof of theorem [7} From lemma El we find, for k > 1 

w v 5 ft(2A; - 1) II j= (N -g-2j)(N-h- 2j) 



2fc-l 



(22) = ~ ^ TT i 

V ; /VCAT — (Oh — 9^ 111 



iV AT. 
0/i(2Jfe -1) ^r 2 / 5-1 \ ( h 



N(N - (2k - 2)) 1JL ^ jy — (2j + 1)/ \ iV - 2j 
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Now fix m, with < 2m < N, and define 
(23) A^A-^Vl-^M X + =U " 



N-mJ\ N-mJ \ N J\ N 

For all k such that 1 < 2k — 2 < 2m, we have the inequalities 

k— 2 

(24) A*- 1 < TT (l 9 ~ l \ ( 1 - — — ) < A 

1 ' " }iV N-(2j + l)J V N-2j! 

which, together with (|22j) . give 

(25) " A - < < P 2fc-l> < ^ _ m) A 

Summing over fc, we obtain the bounds 



k-1 



(26) ^_ < < ^ 



k=l 

where 



fe=i v y fc=i 

Let now m = m(N) be a positive integer sequence such that 

lim m(N) = oo m = o(N). 

N— >oo 

Then A± ~ A and applying lemma [21 we obtain, as iV — ► oo 

(27) ~ K + ~ |^S m ~ [1 - e -™( mK + 1)] . 

We now fix a non-negative real number x. For the convergence of the expression above, we 
specialize the sequence m{N) as follows 

(28) mn = x or m{N) = [xz(N)/2\ 

where n and z were defined in equations (j2Tj) and (fl~2l) . respectively. Due to ([5]), m(N) = o(N), 
and for all sufficiently large iV, we have < m(N) < N/2. Hence (I26D is valid, and our result 
now follows from ([20j) and (J2ZD, noting that (xz(iV) + l)/2 ~ xz(N)/2. □ 

The sequence m{N) given in (|28|) could be defined by the more general requirement that the 
sequence m{N)K(N) converge to a positive real number c, where the limiting cases c = 0, oo 
are excluded as they would lead to trivial distribution functions. We have set c = 1, for 
normalization. We remark that the choice of this constant does not affect the coefficient of 
the distribution function, which determines its limiting value for large arguments. This fact 
will become relevant in the proof of theorem [T] below. 

We now turn to even symmetric cycles, of period t = 2k. Each cycle of this type corresponds 
to an orbit originating on Fix G (or Fix H) , and returning to it for the first time after k 
iterations of L. The arc associated to such a cycle (see equation (|14p ) is a (k + l)-arc, with 
the property that x\ and x^ + i are in FixG. Therefore x\ = y\ and = 2/fc+ij and the 
total number of arcs of this type is equal to N—/2, where the denominator 2 prevents arcs 
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beginning and ending on FixG from being counted twice. The denominator 2 cancels with 
the 2 in the period 2k, so we obtain 

(29) (P£, g ) = ^N^£ 2k k>l 

where 

#E(g-2,h,N-2k) 



^ #V(g,h,N) 
( g(g-l)(N-h) 
N(N-lf 

k—2 k 1 

9 -^^J\(N-g-2j)\{{N-h-2j) iffc>l. 

3=0 j=0 



(30) 



Putting together equation (f29j) and (J30J) , we obtain, for k > 2 



/ P w \ = g(g ~ ^ rf (i 9 — ^ TT (1 - 

\ 2k,gt N ( N _ ( 2fc _ ^) JUL ^ N _ ( 2 j + I) J JUL ^ TV - 

to be compared with equation f)22[) . 

Exchanging g by /i in the formulae above gives the twin quantity P^uu for 2/c-cycles with 

(s e) 

two points of FixH. The distribution function 7£jy' for even symmetric cycles is obtained 
by adding the contributions from the two fixed sets 

[xz/2] 

(31) ^\x) = £ ((P$ 9 ) + (Ht)) 

k=0 

An analysis very similar to that of theorem |4] leads to the following result, whose proof we 
omit for the sake of brevity. 

Theorem 6. For any sequences g, h satisfying f5|) ; with z as in Mty) . we have, as N — > oo 

,(»,e)^ 9 2 + h 2 
(9 + hf 



Proof of theorem [TJ From theorems [Hand [6] we deduce that the probability distribution 
associated to symmetric cycles is given by 

(9 + hf 

Because 

lim 1Z{x) = 1 



1l%' 0) (x) + K^ e) (x) ~ ,7,^ ft(x) = ft(s) . 



this distribution function is properly normalized (cf. remarks following the proof of theorem 
E|), and so it accounts for all cycles. □ 

It is implicit from theorem [T] that, asymptotically, asymmetric cycles have zero probability. 
Asymmetric i-cycles of L = H o G come in pairs, mapped to one another under G. They 
correspond to i-arcs of the type 

(32) (xi,yi,...,x t ,y t ) y j = G{x j )^x j , x j+1 = H(yj) ^ y jt j = l,...,t 
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i.e., they have It distinct entries. We can now determine the decay rate of the probability of 
asymmetric cycles. 

Theorem 7. Consider a composition of random involutions on a set of N points, whose fixed 
sets satisfy As N — ► oo, the measure of the set of points that belong to asymmetric cycles 
is asymptotic to (g(N) + h(N))~ 1 ; in particular, almost all cycles are symmetric. 

Proof. Recall from (I13p that pj a ' is the probability of finding an asymmetric t-cycle. We 
compute its average over E: 

,„(«)v N^ #E(g,h,N-2t) 
{ * ' N #E(g,h,N) 



t-i 



N f£>\ JV-(2j + l)A N-2j 



From the above discussion, the right-hand side of this equation involves counting t-arcs like 
([32]) . Note that the period of the orbit is not present as a factor. This is because it is absorbed 
by the possibility to create i-arcs that are one and the same from (|32p by cyclic permutations 
of the {xj, yj} pairs (t possibilities) and the switch Xj «-> yj together with reversing the order 
of the switched pairs (2 possibilities). 

Let be the associated distribution function 

= E<i> t (B) >. 
t=i 

We now proceed as in the proof of theorem, with A± and z given by (|23p and (|12p . respectively. 
Summing the expression above over t in the appropriate range gives the bounds 

K- < TZ^\x) < K+, 

where 

[xz(N)i [xz(N)i 

k_4 v a 1 n + = 4- y A* 

JV ^ + N ^ + 

t=i t=i 

This time the bounds are geometric progressions; as ./V — > oo, we find 

l-e~ x 



ff (iV) + Miv) ' 

as desired. □ 

The theorems of this section have some immediate consequences. Consider the sequence 
z(N), defined in (|12p . Its denominator g(N) + h(N) is twice the number of symmetric 
cycles [20, proposition 1], and from theorem [TJ asymptotically, the symmetric cycles have full 
measure. These considerations give the following characterization of z(N). 

Corollary 8. Asymptotically, the sequence z(N) defined in $12\) represents the average cycle 
length of the map H o G. 

An immediate consequence of theorems U] and [6] is the following. 

Corollary 9. Let g,h be as above. If g(N) = o(h(N)) or h(N) = o(g(N)), then almost all 
cycles of H o G have even period. 
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3. Occurrence of repeated periods 

This section first gives the proof of theorem B, stated in the introduction. The asymptotic 
relations © will be assumed to hold, together with the condition h = 0(g). The latter 
does not entail loss of generality, as we may always exchange g and h, which amounts to 
considering the inverse mapping L^ 1 . After the proof, we compare the results of theorem B 
with the mechanism that determines the repetition statistics for small period in the Henon 
map ©. 

Proof of theorem B. Let i) be the probability that a randomly chosen (G,H) £ E 
has exactly i cycles of period t. We wish to determine the asymptotic behaviour of \x\ due to 
theorem [7j it will suffice to consider symmetric cycles only. 

We introduce some notation. Let 7 be a i-cycle in f2, and let T be any set of disjoint 
t-cycles. We stipulate that the first point of each cycle belongs to some symmetry line; then 
there are 



ways of choosing such cycles. We define 

(35) £ 7 = {(G, H) G E : 7 is a t-cycle of (G, H)} E r = f] E 1 . 

Note that if the cycles of T were not disjoint, then would be empty. 

Given T, with #T = i, we are interested in the number of pairs in E that support the 
cycles of T, and that do not have any t-cycle in the complement f2 \ T. The parameters of 
the involutions acting on $7 \ T are determined by the type of i-cycles being considered (odd 
period symmetric, even period symmetric on FixG, etc.). Accordingly, we'll consider the 
following sets of pairs of involutions 



E(g — i, h — i, N — it) symmetric, odd period 

(36) E(g — 2i, h,N — it) symmetric, even period, on Fix G 
E(g, h — 2i, N — it) symmetric, even period, on Fix H 

whose cardinality is given by lemma[2j Whenever we don't require specialization, we'll denote 
any of the above sets by the common symbol Ej. The pairs of involutions in Ej act on the 
space Slj with N — it elements. 

Let Ej be the set of elements of Ej that have no i-cycles at all. By the inclusion-exclusion 
principle, we find 

#E, = # I Ei \ |J E 1 I = ^(-l)# r #£ r 
\ ~/en J r 

(37) = £ #E r . 

n>0 #r=n 

The number of non-zero summands is clearly finite. 

To compute the inner sum for fixed t and i, we first choose n distinct cycles, and this can 
be done in (N — it)—/n\ ways. Then we let the involutions act freely on the remaining space, 
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which gives #E i+n possibilities. Thus 



^ #Er = — - #E,+„. 



#r=n 



n! 



Considering the equations (f34"|) and ([37]) , we obtain 



(38) 



fj,(t,i) 



#Ej _ ^ (JV-jt)g#Ej 
^F, n 7:1 ^ l J n ! #E 



i! #E i! 



n=0 



where M = [iV/t - i\ . 

We specialize formula (|38p to the case of odd period i, using the notation For this 

purpose, we consider the first expression E in ([36]) . to obtain 

t) = (Ar " #E(5 " ( * + n) > h ~ ( * + n)) N ~ (i + n)t) 



n=0 



nl 



#E(g,h,N) 



Inserting the values of #E from lemma [21 we obtain, after some manipulation 



M 



(39) 
where 

(40) A {o) (i, 



i\ ^-^ nl 

n=0 



n) 



..i+n ui+n 



n) 



(41) B {o \i,n) 



(N - (i + n)(t - l)pZ 
1 

(i+n)(t-l)/2-l 



n 

3=0 



5-1 



N - (2j + I) 



N-2j 



if t = 1 
if t > 1. 



With a similar argument we derive the probability // e ) for even period 



(42) 

where 
(43) 

(44) 



(e) 



Z z — ' TV 

n=0 



-4 (e) (i,n) 
£ (e) (i,n) 



2(i+n) , h 2(i+n) 



1 + h 



5-1 



2 (JV — (i + n){t — 1))— 

(i+n)(t-2)/2-l , 

(i+n)t/2-l 

x n " 



iV-2j 



/?. 



iV-7 

j=(i+n)(t-2)/2 x 

As in the analysis of section [21 this expression takes into account cycles on both Fix G 
and Fix H, while the factor of 2 compensates for the double-counting that results from the 
presence of two points on each symmetry line. 
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We now consider the asymptotic behaviour of n(t,i), for fixed i. Let y be given by ([7]), 
which we rewrite as 



(45) 



y = x - ln(/) 



x(t) 



(t-l)(g + h) 
2N 



For the purpose of scaling, the period t will be allowed to vary with N. 
Let us consider the expression 



(46) U(a,g,h) = 1 



9 



N -a 



J/2 



1 



h 



N-a 



J/2 



J = (i + n)(t - 1) 0<a<N, 



which, for a = o(N), provides the dominant contribution to the products B(i,n), in both 
cases. A straightforward calculation gives 



(47) n(a,g,h) 
In view of the bounds 



-(i+n)x 



IT(J, g, h) < B(i, n) < 11(0, g-l,h) 



we compute 



U(0,g-l,h) 
U(J,g,h) 



-(i+n)x 



-(i+n)x 



l + O 



n 2 x 2 \ /nxg\ 



+ 



V N J 



We obtain 

and hence 
(48) 



\B(i,n) - e - {i+n)x \ < |n(0 j5 - l,h) -IL(J,g,h)\ 



B(i, n) = e 



-(i+n)x 



l+0{ m )+0 (^ 



This estimate is valid for both odd and even period. 

Next we consider the products A(i,n) in (|40p and (|43p . Recalling the definition of f(N) 
given in equation ©, and the fact that h = 0(g), we have the following asymptotic behaviour 



(49) A(0,0) = 1 A(i,n) = f(N) i+n 

We must consider two cases. 



l + O 



z + n > 



Case J: f(N) -/-> 0. From (|48p and (|49p . we obtain, for both even and odd period 



z! ' — ' n! 

n>0 



n x 



1 + 0(^) + 



n 2 x 2 



9 



(50) 



(e-yy 



( e ~yy 



ra>0 



-l) n e- 



-9V 



+ 



n! 



E 



(-i) 



i! ^— ' n! 

n>0 



+ 0[ - ) + 

.9, 



14 JOHN A. G. ROBERTS AND FRANCO VIVALDI 

where we have used the identities 

V<_i£_=0 = O<l) Ttll2L = - 1 - = o m . 

^ n\ ^ n\ e 

n>0 n>0 

Now we adjust t so that y is constant. This gives x = 0(ln(/)), and for both odd and even 
period the error terms in (|50p vanish. We have convergence to a Poisson distribution 



(51) n(t,i) ~ e~ a -r a = e-^. 

i! 

Note that for the above result to hold, we have only required that f(N) -/-> 0, not the existence 
of a limit for f(N). 

If /(iV) tends to a non-zero finite limit c, there is no logarithmic term in the function y(N) 
— cf. (|45p — and it is possible to consider the probability fj,(t, i) without scaling the period t. 
For fixed t, we have x = O(gfN) — ► and y — > — ln(c), and equation ([50]) gives 

& 

H(t,i) ~ e c — , 
z! 

independent of t. 

Case II: f(N) -» 0. From ggj we have , ti) ~ ^i+n (Kronecker's delta), and hence 

The proof of theorem B is complete. D 

We remark that theorem A can actually be derived from theorem B since the distribution 
TZ(x) follows from first principles as 

[xz\ [N/tj 

(52) K(x) = lim lZ N {x) = lim V - V Oi 

N— >oo N^oo — » iv ■ £ — ' 

t=l i=l 

equivalently (Pt) follows from 

(53) < P *> = ^ E 

i=i 

Note that the variable x in ([35]) is, asymptotically, the variable x in IZ(x). Taking the result 
([ST]) gives 

oo oo c 

cr 



J^i //(*,*) = ae 



00 (i-i) 



i=i i=i v 1 

Since e~ y = f e~ x from (|45|) . we have 

(54) <P t )~l/e-*. 

If t = 2/c — 1 is odd, we substitute for / from ([6]) and x from ([4"5]) to see 

(55) (P 2fc -i)~^|^ e -^. 

Comparing this to ([25]) and noting there that A± _1 ~ A fc_1 = e^ -1 ) 111 ^ ~ e -(k-i)n s \^ ows 
that we will recover the result of theorem [3] A similar exercise on (|54p for t = 2k using the 
appropriate / from ([6]) yields the result of theorem [6] and hence the distribution TZ(x) upon 
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summation of odd and even contributions. We have chosen to give the independent proof of 
Theorem A in section [2] as it is simple and direct and bypasses the calculation of n(t,i). 

The applicability of Theorem B to algebraic maps over finite fields should be considered 
with care. For fixed t, the maximum number of i-cycles for an algebraic map is plainly finite, 
because the points in a cycle are solutions of a set of polynomial equations; therefore, in the 
large field limit, one cannot possibly expect convergence to ([8]). Under what circumstances is 
then theorem B a model for such maps? 

To articulate this question, let us consider a specific example, namely the Henon map ([2]), 
with a = 1. The set Fix G is the line y = x, and so over the field ¥ p , this map has parameters 
g = h = p, hence a = c = 1 in theorem B. Let us consider symmetric cycles, which, due to 
theorem A, are expected to dominate the dynamics. When the fixed sets are one-dimensional, 
the study of symmetric orbits is greatly simplified, because the points in these orbits are roots 
of polynomial equations in one indeterminate (as opposed to being points on varieties). Now, 
any symmetric odd cycle has a point (x, x) on Fix G, and the values of x for the period t = 5 
are found to be the roots of the irreducible polynomial 



Thus the number of symmetric 5-cycles over the field ¥ p is equal to the number of linear 
factors in the factorization of modulo p. In particular, the maximum number of 5-cycles 
is equal to 6, the degree of $5. 

In general, polynomial factorization over a finite field is a highly irregular process, subject 
only to a slight deterministic constraint — Stickelberger's theorem — which determines the 
parity of the number of factors [21 section 4.8]. Significantly, this process is ruled by a prob- 
abilistic law — Cebotarev's density theorem — which determines the probability associated 
with each factorization type (number and degree of the factors) in terms of properties of the 
Galois group of the polynomial [15\ p. 129]. Specifically, for each factorization with distinct 
roots, the degree of the factors defines an additive partition of the degree of the polynomial. 
Each partition in turn identifies uniquely a conjugacy class of the Galois group, with the 
terms in the partition corresponding to cycle lengths in the cycle decompositions of the per- 
mutations in the class. The probability of a factorization type is then given by the relative 
size of the corresponding class in the Galois group. If the Galois group of $ is the symmetric 
group, the largest possible one, then every factorization type does actually occur, and with 
positive density among all primes. 

It can be verified that the Galois group of $5 is in fact Sq, the set of all permutations 
of deg(<]?5) = 6 objects [3]. There are therefore 11 possible types of factorizations of <I>5, 
corresponding to the 11 additive partitions of the integer 6. We factor ^(x) modulo a few 
consecutive primes 



These factorizations correspond to the partitions 1 + 1 + 1 + 1 + 1 + 1 = 6 = 1 + 1 + 4, 
respectively. We see that, over the field ¥ p , the Henon map has six symmetric 5-cycles for 
p = 6563 (the smallest prime for which this happens), no 5-cycles at all for p = 6569, and 
two 5-cycles for p = 6571. 

The degree of the polynomials &t grows exponentially with t, and it is reasonable to expect 
that, typically, they will be irreducible, with large Galois groups. For the symmetric group, 



$ 5 {x) 



x 6 - 2 x 5 + 5 x 4 - 6 x 3 + 8 x 2 - 4 x + 3. 



$ 5 (x) 



' (x + 508) (x + 917) (x + 1165) (x + 1486) (x + 3168) (x + 5880) 
x 6 + 6567 x 5 + 5 x 4 + 6563 x 3 + 8 x 2 + 6565 x + 3 
k (x + 1265) (x + 3889) (x 4 + 1415 x 3 + 2999 x 2 + 4200 x + 2683) 



(mod 6563) 
(mod 6569) 
(mod 6571). 
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the Cebotarev's probability u(t, i) of occurrence of % symmetric cycles of period t is just the 
probability of having exactly i fixed points in the appropriate symmetric group, namely 

(56) ^ t ,i) = ^-J2(-iy- d = deg* t (x). 

*' 3=0 J - 

In the large t limit, v(t,i) converges to e _1 /i!, independent of t. This agrees with the predic- 
tion of the random involution model. 

It is instructive to compare the data (|56|) from Cebotarev's theorem with the prediction of 
theorem B, for the period t = 5. 



i 


Cebotarev thm 


thm B 





53/144 = 0.3681... 


0.3679. 


1 


11/30 = 0.3667... 


0.3679. 


2 


3/16 = 0.1875... 


0.1839. 


3 


1/18 = 0.0555... 


0.0613. 


4 


1/48 = 0.0208... 


0.0153. 


5 





0.0031. 


6 


1/720 = 0.0014... 


0.0005. 



We see that, even for such a low period, for i < 3 the agreement is already reasonable. 

Thus, in the two-dimensional case, some statements on the probability of multiple occur- 
rences of symmetric i-cycles translate into statements on the maximality of the Galois groups 
of &t- Even though this formulation would allow one to re-cast some dynamical questions in 
algebraic terms, we note that problems related to Galois groups of iterated polynomials are 
often quite difficult. 

4. Concluding remarks 

When the reduction over a finite field of a reversible algebraic map with a single family of 
reversing symmetries is a permutation, numerical experiments suggest the period distribution 
function is given by the period distribution function lZ{x) of ((4]). In this paper, we have shown 
that lZ{x) corresponds to the expected period distribution function in a probability space of 
pairs of involutions characterized by the cardinalities of their fixed sets. We find that the 
same probability space furnishes a Poisson law for the probability of cycles with the same 
period, also in agreement with numerics. 

We conclude with some remarks: 

(1) To prove theorem A, only mild regularity conditions — the existence of the limits ([5]) - 
were imposed on the behaviour of the sequences g(N), h(N), specifying the cardinalities of 
the fixed sets of the two involutions, as the size N = q n of the phase space increases (with n 
the dimensionality of the map). Similarly, theorem B involves checking the behaviour of the 
sequence f{N) in (|6|). We point out that when a reversible map of a finite set is obtained by 
reducing an algebraic mapping to a finite field, the behaviour of such sequences is strongly 
constrained, resulting in algebraic growth. 

A simple but significant case comprises maps for which the fixed set Fix G is an affine 
subspace of some linear space. An example is given by two-dimensional reversible polynomial 
automorphisms over some field K, where the existence of a normal form for the two involutions 
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ensures that their fixed set is a single point or a line pQ. Localization to the finite field ¥ q 
gives a phase space with N = q 2 points, and since a line in the affine plane ¥ 2 has q points, 
we have an algebraic growth: g(N) + h(N) = 2q = . In dimension d, each of FixG and 
Fix H can have integer dimension ranging from to d — 1, and if these sets are affine spaces, 
we obtain the sequence g(N) + h(N) = N r + N s , with r, s € Q, < r, s < 1. The conditions 
([5]) are then met as long as one of r or s is nonzero, i.e. as long as one of Fix G and Fix H is 
at least one-dimensional. Furthermore, in this situation, we see from ([6]) that f(N) = N r+S ~ 1 
or f(N) = (N 2 ^ 1 + N 2f >- 1 )/2. So either / = 1 (e.g., when r = s = 1/2 ^g = h = ^/N) oi 
f — > oo algebraically or / — > algebraically. 

More generally, the sets Fix G and Fix H are algebraic varieties over ¥ q determined by 
n polynomial equations, and we are interested in determining the number of points on these 
varieties. The number M of points on an r-dimensional irreducible algebraic variety in the 
n-dimensional projective space of degree d over the field ¥ q satisfies the Hasse-Weil bound 

m 

(57) \M - q r \ <(d-l)(d- 2)q r ~ 1 / 2 + Cq r ~ x 

where C = C(n,r,d). For r = 1 we obtain a square root bound, and we can also replace 
(d — l)(d — 2) by the genus of the curve. The result (|57p ensures that, in the large field limit, 
there is algebraic growth of the cardinality of fixed sets. So the discussion of the previous 
paragraph applies concerning the validity of the asymptotic conditions © and the behaviour 
of f(N). 

(2) Experiments show that algebraic maps without time-reversal symmetry behave quite 
differently from reversible ones, when represented on a finite field |20j. In the absence of 
reversibility, the cycle distribution is found to be that of a random permutation, given by the 
identity function on the unit interval. This distribution, obtained via the scaling z(N) = N, 
has led to a general conjecture [201 conjecture 2]. 

(3) When a general rational reversible map is reduced over finite fields, we expect that points 
where the denominators vanish (singularities) will occur. Restricted to its periodic points, 
the reduced map is still a permutation. Numerically, we find the distribution of periods is still 
governed by 7Z(x), consistent with Theorem A [9l[22]. But nonperiodic orbits will also appear 
in the reduction of such a map. We are presently investigating the extension of the random 
involution model to explain various aspects of the dynamics of the reduction to finite fields 
of general rational reversible maps [22]. We are also extending our numerical experiments to 
cover the asymptotic regime ¥ q , q = p k , k — ► oo. 
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